Search for: All records

Abstract Mineral dissolution releases ions into fluids and alters pore structures, affecting geochemistry and subsurface fluid flow. Thus, mineral dissolution plays a crucial role in many subsurface processes and applications. Pore‐scale fluid flow often controls mineral dissolution by controlling concentration gradients at fluid‐solid interfaces. In particular, recent studies have shown that fluid inertia can significantly affect reactive transport in porous and fractured media by inducing unique flow structures such as recirculating flows. However, the effects of pore‐scale flow and fluid inertia on mineral dissolution remain largely unknown. To address this knowledge gap, we combined visual laboratory experiments and micro‐continuum pore‐scale reactive transport modeling to investigate the effects of pore‐scale flow and fluid inertia on mineral dissolution dynamics. Through flow topology analysis, we identified unique patterns of 2D and 3D recirculating flows and their distinctive effects on dissolution. The simulation results revealed that 3D flow topology and fluid inertia dramatically alter the spatiotemporal dynamics of mineral dissolution. Furthermore, we found that the 3D flow topology fundamentally changes the upscaled relationship between porosity and reactive surface area compared to a conventional relationship, which is commonly used in continuum‐scale modeling. These findings highlight the critical role of 3D flow and fluid inertia in modeling mineral dissolution across scales, from the pore scale to the Darcy scale. 

Abstract Fluids with different densities often coexist in subsurface fractures and lead to variable‐density flows that control subsurface processes such as seawater intrusion, contaminant transport, and geologic carbon sequestration. In nature, fractures have dip angles relative to gravity, and density effects are maximized in vertical fractures. However, most studies on flow and transport through fractures are often limited to horizontal fractures. Here, we study the mixing and transport of variable‐density fluids in vertical fractures by combining three‐dimensional (3D) pore‐scale numerical simulations and visual laboratory experiments. Two miscible fluids with different densities are injected through two inlets at the bottom of a fracture and exit from an outlet at the top of the fracture. Laboratory experiments show the emergence of an unstable focused flow path, which we term a “runlet.” We successfully reproduce the unstable runlet using 3D numerical simulations and elucidate the underlying mechanisms triggering the runlet. Dimensionless number analysis shows that the runlet instability arises due to the Rayleigh‐Taylor instability (RTI), and flow topology analysis is applied to identify 3D vortices that are caused by the RTI. Even under laminar flow regimes, fluid inertia is shown to control the runlet instability by affecting the size and movement of vortices. Finally, we confirm the emergence of a runlet in rough‐walled fractures. Since a runlet dramatically affects fluid distribution, residence time, and mixing, the findings in this study have direct implications for the management of groundwater resources and subsurface applications.